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The local electronic structure of Nd2Cu04 is determined from ab-initio cluster calculations in 
the framework of density functional theory. Spin-polarized calculations with different multiplicities 
enable a detailed study of the charge and spin density distributions, using clusters that comprise up 
to 13 copper atoms in the Cu02 plane. Electron doping is simulated by two different approaches and 
the resulting changes in the local charge distribution are studied in detail and compared to the cor- 
responding changes in hole doped La2Cu04. The electric field gradient (EFG) at the copper nucleus 
is investigated in detail and good agreement is found with experimental values. In particular the 
drastic reduction of the main component of the EFG in the electron-doped material with respect to 
La2Cu04 is explained by a reduction of the occupancy of the Zd^,,-i_^2 atomic orbital. Furthermore, 
the chemical shieldings at the copper nucleus are determined and are compared to results obtained 
from NMR measurements. The magnetic hyperfine coupling constants are determined from the spin 
density distribution. 

PACS numbers: 74.25.Jb, 74.62.Dh, 71.15.Mb, 74.72. Jt 



I. INTRODUCTION 

Soon after high-Tc superconductivity was dis- 
coveredii^ in hole-doped La2_a;(Ba,Sr)a;Cu04 and 
YBa2Cu307_i- materials, it was found^ that electron- 
doped Nd2-xCca:Cu04 also supcrconducts. All these 
substances contain Cu02 planes separated by block lay- 
ers containing rare-earth ions, and, moreover, the par- 
ent compounds are insulating and exhibit quasi two- 
dimensional antiferromagnetism. The generic phase dia- 
grams of electron- and hole-doped cuprates, however, are 
somewhat different and are not totally symmetric with 
respect to hole or electron doping. This is not unexpected 
since there are quite significant structural differences, fn 
La2Cu04, the copper is six- fold coordinated with oxygens 
whereas in Nd2Cu04, the apical oxygens are missing. To 
gain insight into the mcchanism(s) of high-Tc supercon- 
ductivity, investigations on the origin of similarities and 
differences between hole- and electron-doped materials 
are very important. 

In a simplistic ionic picture the valencies of the Cu 
and the O arc the same in the insulating state of the 
parent compounds La2Cu04 and Nd2Cu04. The former 
substance becomes metallic upon hole doping and it is 
assumed that the doped hole goes to the planar oxygen 
sites. On the other hand, in Nd2Cu04, the metallic state 
is achieved by electron doping and the doped electrons 
are assumed to reside on the copper sites. 

Electronic structure calculations^ in the framework 
of the full-potential linearized augmented-plane-wave 
method within the local density approximation have 
shown very similar band structures for Nd2Cu04 and 
La2Cu04. Details of the influence of electron or hole 
doping in these materials, however, cannot be obtained 
from band structure calculations except in the approxi- 
mation of rigid bands. In contrast, cluster calculations 
are ideally suited for investigations of how the local elec- 



tronic structure is influenced by doping. Using two dif- 
ferent methods of simulating doping, we have previously 
showni that simulations of hole doping in La2Cu04 pro- 
vide insight into the changes in the population of various 
orbitals on the copper and oxygen atoms. Both methods 
yielded on the whole the same results and the changes 
calculated for the electric field gradient (EFG) at the 
copper agreed with experiments. Parts of this work deal 
with an adaptation of these methods to the case of the 
electron-doped cupratc Nd2Cu04. 

A large body of data obtained by nuclear quadrupole 
resonance (NQR) and nuclear magnetic resonance 
(NMR) gives detailed information on the changes 
of the EFG values in all high Tc superconduc- 
tors. Several NMR experiments were performed in 
the electron-doped Nd2Cu04 compoundSi2i2i2iiSiii, the 
Pri_2:LaCea;Cu04_5 materiali^, the infinite layer com- 
poun d-^^d'^d^ Sri_a;La2;Cu02, and its parent compoundi^ 
SrCu02. The first surprising feature of these results is a 
relatively low ^'^Cu NQR frequency {vq) measured in the 
undoped parent compound compared to the hole doped 
compounds. In the undoped parent compounds the 
quadrupolar frequency has been measured to be 14 MHz 
for the the Nd2Cu04 compound and 7.4 MHz in the infi- 
nite layer. This is less by more than a factor of two than 
in La2CuO4i33.0 MHz (Ref.IIl)) and YBazCugOy (31.5 
MHz (Ref. n^) ) compounds. The second striking feature 
of all these experiments is the extreme sensitivity of vq 
on the doping level. vq\s<2 MHz for Ndi .8.5Cen.i.qC u04 
(Ref. 113), < 0.5 MHz in Pro.9iLaCeo.o9Cu04-5_(Ref.|l3) 
and finally < 3 MHz in Sro.9Lao.iCu02 (Ref.[lJ). This 
strong doping dependence of the field gradient in the 
electron-doped compounds has been explained by a sim- 
plistic ionic model. It was shown to be due to fully occu- 
pied 3dj.2_y2 orbitals which give no contribution to the 
field gradient if all 3d orbitals are occupied. 

In the present work we report on cluster calculations of 
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the local electronic structure in Nd2Cu04 and its changes 
upon electron doping. The results are compared to those 
for hole-doped La2Cu04. The contributions to the EFG 
are studied in detail and it is explained why the EFG 
values in Nd2Cu04 are much smaller than in the hole- 
doped substances. 

In Sec. we present the details of the technique of 
the cluster calculations. In Sec. IIIII the results for the 
electronic structure including the charge distribution, 
the EFG and their comparison to NMR experiments are 
given. Especially, in Sec. IIII El results for the chemi- 
cal shieldings at the copper sites in undoped Nd2Cu04 
are presented and are compared to results obtained from 
NMR experiments. The results for the spin density dis- 
tribution and the calculation of magnetic hypcrfine cou- 
plings are discussed in Sec. IIVI Sec. contains a sum- 
mary and conclusions. 



II. THE CLUSTER METHOD AND 
COMPUTATIONAL DETAILS 

A. General Remarks 

The idea of the cluster method is to solve the many- 
body Schrodinger equation for a portion of the crystal 
which we call the cluster. A cluster is a careful selection 
of a contiguous group of ions within a solid. The spe- 
cific choice of the atoms that make up a cluster is such 
that it allows predominantly localized properties of a tar- 
get atom and its vicinity to be calculated. A property is 
called local if it can be determined by the electronic struc- 
ture of a few neighboring atoms in the crystal. Electric 
field gradients (EFG) and hyperfine fields are such lo- 
cal properties and have been successfully calculated by 
the cluster method (see for example Refsi^Si^S). A clus- 
ter consists of three regions. The target atom and at 
least its nearest neighboring atoms form the center of 
the cluster and the corresponding electrons are treated 
most accurately using first-principles all-electron meth- 
ods. This core region is embedded in a large cloud of 
a few thousand point charges at the respective lattice 
sites imitating the Madelung potential. Point charges at 
the border of the core region are replaced by basis-free 
pseudopotentials to improve the boundary conditions for 
the electrons in the cluster core. These pseudopotentials 
make up the so-called screening region. 

It would be desirable to have clusters that contain as 
many atoms as possible in the core region but there 
are two computational limitations to the cluster size: 
the available computer resources and the convergence of 
the self-consistent field procedure. However, these limits 
have been pushed further since our first use of the cluster 
technique (see Ref. l2l|) which now allows the computation 
of larger clusters and hence improve the bulk-to-surfacc 
ratio. 

In this work three clusters of different size (with 5, 
9; and 13 copper atoms in the core region) are studied 



Cluster NEB 
Cu5 0i6/Cu8Nd24 2 1 29 5 403 
Cu9024/Cui2Nd32 33 483 663 
Cui3 036/Cui2Nd48 49 711 975 

TABLE I: Compilation of the clusters used for Nd2Cu04 and 
some of their defining properties: number of atoms simulated 
by a full basis set (A'^), number of electrons (E) and number 
of basis functions (B). 

and their constitutive properties are given in Table ^ 
For clarity these clusters are labeled X/Y where X is the 
chemical formula of the core region and Y the formula of 
the screening region. 



B. Computational Details 

The atomic positions were chosen according to crys- 
tallographic structure determinations. The crystal struc- 
ture of Nd2Cu04 is the so called T' structure. A charac- 
teristic of this structure is the absence of apex oxygens 
in contrast to hole doped materials. In this case the core 
region consists only of planar Cu and O. The lattice pa- 
rameters (a = 6 = 3.95 A and c = 12.07 A) and the posi- 
tions of the atoms within the unit cell have been adopted 
from Refsi^SiSi. The core region of the smallest cluster 
(CU5) for Nd2Cu04 incorporates five planar Cu atoms 
and their 16 nearest planar oxygen atoms. The screen- 
ing shell includes bare pseudopotentials for 8 Cu^+ and 
24 Nd'^^ ions. The Madelung shell contains more than 
6000 point charges. According to the attributed valences 
of copper and oxygen of 2+ and 2~, respectively, this 
cluster is given a total charge of 22", leading to a to- 
tal number of 295 electrons. We will also use Cug and 
Cui3 clusters whose properties are given in Table |l| The 
largest cluster (Cuia) is schematically depicted in Fig.^ 

For the atoms in the core region the standard 6-31 IG 
basis sets^S were employed. The electronic structures in 
all the clusters described above were determined with the 
Gaussian03 quantum chemistry package^^ in the frame- 
work of density functional theory incorporating the ex- 
change functional proposed by BeckeSSiSi, together with 
the correlation functional of Lee, Yang, and Parr— (spec- 
ified by the BLYP keyword in the Gaussian03 program). 
We note that band structure calculations with the local 
density approximation (LDA) give no magnetic solution 
for Nd2Cu04. In cluster calculations, however, the spin 
state can be chosen and in all cases (irrespective of the 
functional, LDA or generalized gradient approximation) 
the state with the lowest variational energy is found for 
the antiferromagnetic arrangement^Siffi. 

The results of each calculation were examined with the 
MuUiken population analysis which gives a description of 
the charge and spin densities in terms of the individual 
atoms and also the constituent orbitals. Other properties 
such as EFGs and hyperfine fields were also recorded at 
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each copper center. It is important to realize that prop- 
erties calculated with the cluster method as e.g. EFGs 
or hyperfine fields depend only marginally on the cluster 
size (see a detailed analysis in Ref. Isij) . 




FIG. 1: (color online). Schematic representation of the 
Cui3 036/Nd48Cui2 cluster. The 13 small black spheres de- 
note the copper atoms simulated with basis sets, the 12 small 
grey (brown) spheres at the border of the cluster are the cop- 
pers represented by pseudopotentials. 



C. Simulation of doping 

There is an important aspect to consider when study- 
ing the change of the electronic structure upon doping 
cupratcs by means of cluster calculations. If the doping 
is produced by impurity atoms (as for example replacing 
La3+ by Sr2+ in La2Cu04 or by Ce'*+ in Nd2Cu04) 

one can investigate the local changes in the immediate 
neighborhood of the impurity atoms. In this situation the 
doped charge carriers are "pinned" around the dopant 
and are thus interpreted as localized holes (electrons). In 
contrast, we would like to investigate holes (electrons) 
distributed over the Cu02 plane. Therefore we abstain 
from changing the chemical composition of our clusters 
but instead we simulate doping cither (I) by simply re- 
moving (or adding) one electron or (II) by creating hole 
(electron) excess at the center of the cluster by an ad- 
ditional electric field. Note that in both methods, the 
nature of the dopant is not specified. 

In this work wc generalize the methods used in Refsi^iS^ 
for the hole doped material (La2Cu04) to the electron- 
doped material (Nd2Cu04). If we simulate doping by 
adding an electron to the cluster (method I) the spin 
state of the cluster changes. The corresponding doping 
level is about 20 % for a Cus cluster, 10 % for a Cug 
cluster, and 8 % for a Cuis cluster. It is therefore not 
possible to simulate an arbitrary doping level with this 
method. 

For the second method wc introduce additional point 
charges at the periphery of the cluster to create an elec- 
tric field across the cluster in order to move the charge 
towards or away form the target atoms in the cluster cen- 
ter. We call this method the peripheral charge method. 
Note that this added system of charges has no physical 
meaning. It is convenient because it can be continuously 
altered so that the charge can be progressively directed 
to or extracted from the target atoms in the center of the 



cluster. This implies that the cluster must be sufficiently 
large so that ions of interest are not near the cluster edge. 
In contrast to the first method the cluster and the system 
of charges always have the same number of electrons and 
the same spin state, but the charge of the cluster can be 
progressively changed in a manner expected by doping. 
Of course this model has only a limited range of applica- 
bility and would ultimately break down for high doping 
levels since no electron is actually removed or added. The 
peripheral charge method serves to continuously interpo- 
late the results already obtained with the first method 
(with remarkable success, see Sec. lIIIi^ . 

III. ELECTRONIC STRUCTURE 

A. Charge distribution 

In cluster calculations, the charge distribution is de- 
termined by the occupied molecular orbitals (MOs). The 
TO*'* MO is represented as a linear combination 

n n riK 

^,n{r^ = E '^™(^^- = E E c^''SA-,fe(r - Rk) 

K=l K=l k=l 

(1) 

of riK atomic basis functions BK,k centered at the nu- 
clear sites K = 1, . . . , n, and the c^'*^ are the MO coeffi- 
cients^. The charge density is then given by J^m 4'm{r)'^ 
and contains both on-site and overlap terms. In the fol- 
lowing we will use the MuUiken population analysis to 
describe the details of the charge distribution. 

In the left part of Fig. |21 we have represented the 
partial MuUiken populations, pc, for the relevant or- 
bitals: 3d2;2_j,2, Sd^z'i-r^, and 4s for Cu and 2pct for 
Op. All other orbitals have populations that differ only 
marginally (< 0.2%) from 2. Results have been ob- 
tained from calculations with clusters of three different 
sizes: Cu50i6/Cu8Nd24 with spin multiplicity M = 4, 
Cug024/Cui2Nd32 with M = 2, and Cui3036/Cui2Nd48 
with M = 6 for the imdopcd case. We find that the total 
energy is lowest for the multiplicity that corresponds to 
an antiferromagnetic spin alignment. We repeated these 
calculations with one additional electron and, again, with 
a spin multiplicity leading to the lowest ground state en- 
ergy {M = 3, 1, and 5 for Cus, Cug, and Cui3, respec- 
tively). This doping simulation changes the total charge 
in the unit cell. We therefore define a quantity p(3) as the 
sum of the MuUiken charge on the central copper atom 
and that on the adjacent oxygen atoms, 

p(3) = p(Cu) + 2p(0p) + 2, (2) 

which approximately represents the fractional excess 
charge in the unit cell. Ideally, in the undoped sub- 
stance p{3) should vanish. Due to open boundaries in 
the cluster method this is not quite exactly the case but 
since all calculated values arc sufficiently close to zero we 
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have rcnormalizcd the values for p(3) such that they are 
indeed zero in the undoped compounds. The full large 
(red) symbols are the partial Mulliken populations calcu- 
lated in the Cu,5 cluster in the undoped case (p(3) = 0) 
and with one additional electron at p{3) « —0.1. The 
full small (blue) symbols denote pc values obtained from 
a cluster with 13 Cu again in the undoped case and with 
an additional electron at p(3) w —0.05. The results of the 
calculation for the cluster with 9 atoms are not shown al- 
though in the undoped case they are almost identical to 
those of Cui3. However, when adding an electron, the re- 
sults exhibit a peculiar superstructure for which presently 
we do not have a unique explanation. The open symbols, 
finally, represent partial populations that have been de- 
termined with the peripheral charge method described 
in Sec. Ill Ci that allows to simulate the doping with frac- 
tional charges. These values have been obtained with the 
Cu5 (Cuia) cluster and result from calculations that all 
used the same spin multiplicity M = 4 (M = 6) as in 
the undoped case. (Note that the partial populations of 
the 3d322_r2, 4s, and 2pg. orbitals in the Cus and Cui3 
clusters for p{3) = are so close that the correspond- 
ing markers nearly completely overlap.) It is seen that in 
the undoped case the partial populations only marginally 
depend on the cluster size. As a function of doping the 
Pc values calculated for the two differently sized clusters 
(filled symbols) change monotonically and are reasonably 
well interpolated by those obtained with the peripheral 
charge method (open symbols). 

For a general discussion it is instructive to compare 
these results with those obtained previously by StoU et 
al^ for hole doped La2-2;Sr2;Cu04 which are shown in 
the right part of Fig. |21 Note that in this case p{3) is 
defined by 

p(3) p(Cu) + 2p(0p) + 2p(0a) + 6. (3) 

Again, the full symbols denote partial Mulliken pop- 
ulations of the relevant AO calculated for the undoped 
material {p{3) ~ 0) for the cluster with 5 Cu atoms 
with M = 4 and with one electron less and M ~ 3 
{p{3) w 0.12). The open symbols arc obtained from 
calculations with the peripheral charge method. 

By comparing the left and right sides of Fig. |21 it is 
first observed that there are no drastic differences be- 
tween the electron and hole doped materials although 
the latter have apical oxygens whose 2p2 orbitals arc not 
completely filled. The populations of the Cu 3dx2_y2 or- 
bitals (diamonds) and the O 2po. orbitals (squares) which 
carry most of the spin density slightly decrease from left 
to right. There is of course a discontinuity at p(3) ~ 
which is due to the different structures (the lattice con- 
stant a for Nd2Cu04 is 4.8 % larger than for La2Cu04). 

It is remarkable that the Cu Sd^z'^-r^ AO (triangles 
down) in the electron-doped material is less occupied 
than in La2Cu04. This is compensated by a significantly 
larger population of the Cu 4s AO (circles). In Fig.|31the 
density difference in the 4s and Sd^z^-r^ AOs between 



the Cui3036/Cui2Nd48 and Cui3062/Cui2La74 is shown 
which demonstrates the increased population of the 4s 
at the expense of the 3d3^2_^2 AO. 



B. Sensitivity of EFG values on charge distribution 

A detailed description of the evaluation of EFGs in 
the framework of cluster calculations is found in Ref. 13^ . 
Here we just point out the relevant aspects that are con- 
nected to the rearrangement of charge distributions due 
to electron doping. 

We assume in the following that the target nucleus Kt 
is at Rkt = 0- The contribution of the MO 4>m to the 
EFG at Kt is given by the matrix element 

= { 0™(f)| I cf>^{r) ) 

n n 71 K nL 

= X] X] X! X! "^m '''^m' ^ 
A'=l L=l k=l 1=1 

x{BKAr-RK)\ ^"^^^ y ^' \BL,i{r-RL))- (4) 

This matrix element contains contributions from basis 
functions centered at two nuclear sites K and L. Thus 
we can identify three types of contributions: (i) on-site 
terms from basis functions centered at the target nucleus 
[K = L = Kt, contribution I), (ii) mixed on-site and 
off-site contributions (II), and (iii) purely off-site terms 
with K ^ Kt and L ^ Kt (HI). 

In addition, there is a contribution coming from all 
nuclear point charges Zx which we denote by j . The 
charges of the bare nuclei at sites K ^ Kt are screened 
by the matrix elements with K = L. Therefore the com- 
bined contributions from III and the nuclei (W^) are 
small. The partitioning of the contributions to the EFG 
tensor Vij thus reads 

V,, = 'V,j + "V,, + '"V,, + W,, . (5) 

More details about these regional partitions can be found 
in Ref. 113 

We collect in Table Hll the various contributions to Vzz 
at the copper obtained for the Cui3036/Cui2Nd48 clus- 
ter for both the undoped and the electron-doped calcu- 
lations. In the undoped case, the total value Vzz = 0.538 
is obtained from contributions I (0.428), II (0.160), III 
+ nuclei (—0.015) and a small correction due to the sur- 
rounding point charges (—0.033). In the electron-doped 
case the total EFG of Vzz ~ 0.235 is substantially smaller 
which is mainly due to a diminished contribution from 
region I. Focusing on region I we note that the added 
contributions from the three AO idxy, idzx, and 3dj,2 
are —0.072 (undoped) and —0.064 (doped). 
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FIG. 2: (color online). Partial MuUiken populations of the atomic orbitals versus doping level p(3) for Nd2-iCea;Cu04 (left 
panel) and La2-iSriCu04 (right panel). The following orbitals are represented: 3d3^2_^2 (triangles down), 4s (circles), apical 
oxygen 2pz (triangles up), 3d^2_y2 (diamonds) and the planar oxygen 2p^ (squares). Open symbols refer to peripheral charge 
calculations and filled symbols refer to calculations where electrons are added or removed as appropriate (method I). 



undoped electron-doped 







(M = 6) 


(M = 5) 


I 


P 


-1.431 


-1.433 




d^2_y2 


-6.378 


-6.613 






-9.128 


-9.088 




d3^2_^2 


8.309 


8.222 




dxz 


4.528 


4.512 




dyz 


4.528 


4.512 


II 


S 


0.374 


0.396 




P 


0.013 


0.013 




d^2_y2 


-0.021 


-0.020 






0.001 


0.001 




d3z2_r2 


-0.207 


-0.217 






0.000 


0.000 




dyz 


0.000 


0.000 


III + Nuclei 




-0.015 


-0.017 


Point charges 




-0.033 


-0.033 


Total 




0.538 


0.235 



TABLE II: Contributions I, II, III to the EFG component 
Vzz for the central copper (in atomic units), calculated for 
the Cui3 cluster representative of Nd2Cu04. 



A careful analysis of Table |n] shows that the most sig- 
nificant difference is due to the dx2_y2 and the d^^'^-r^ 
parts of contribution I. Only very small corrections are 
obtained by the other contributions. The influence of 
part III of the electrons on the other atoms is nearly 
fully compensated by the nuclear charges of these atoms 



and contributes only by a very small amount. 



To understand the changes of the EFG on the copper 
sites upon doping, it is therefore sufhcient to concentrate 
on the MOs which contain partially occupied 3d2.2_j,2 
and 3d3^2_y.2 AOs at the target nucleus. The EFG is 
then given by 



(r )^2_j^2), 
(6) 

where R is the contribution from all other orbitals and 
arc the partial occupation numbers of the corresponding 
orbitals which are gathered in Table ITTIl 

The fact that the partial occupation numbers A^ of 
the AOs are very similar to the partial MuUiken popu- 
lations pc which are also included in Table IIIII gives evi- 
dence that despite their simplistic concept the MuUiken 
charges can be significant quantities in the analysis of 
charge distributions in cuprates. In particular, the dif- 
ferences N3z2_^2 — Nj.2_y2and Pci'id^z^-r'^) —Pc{-idx2_y2) 
are very similar. In Sec. lIII Pl we will therefore use a sim- 
plified version of Eq. J^Jl for a more general analysis of 
the Cu EFGs in cuprates. 
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If we adopt the value of the quadrupole moment ^^Q = 
—0.211 b derived by Sternheimer— we obtain from our 
calculations a quadrupolar frequency oi^^iyg = 13.3 MHz 
in the undoped Nd2Cu04 which is in very good agree- 
ment with the value of i/q = 14 MHz reported by Abe 
et aim'. In the electron-doped materials we get ^^i^q ~ 
5.8 MHz in the Cui3 cluster (corresponding to a doping 
level of about 8 %) and ^"^i^q = 3.6 MHz in the Cu5 clus- 
ter (corresponding to 20 % doping) . Experimental values 
for quadrupole frequencies in doped materials are gener- 
ally much smaller than in the undoped case and close to 
zero {vQ < 2 MHz in Ref. ITol for a doping level of 15 %). 
This reduction compared to the undoped substance is at 
least qualitatively reproduced by our calculations. 



undoped electron-doped 

3dj.2_y2 3dr^22_^2 3dj.2_y2 Sd^^^2_j,2 

N 1.379 1.820 1.435 1.806 

Pc 1.430 1.865 1.483 1.856 

(r~^) 8.094 7.991 8.064 7.968 

TABLE III: Partial occupation numbers N, partial MuUiken 
populations pc, and expectation values (r~'^ ) for the relevant 
copper AOs in undoped (M = 6) and doped (M — 5) Cuis 
clusters representative of Nd2Cu04. 
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FIG. 3: (color online). Electron density difference contour 
plots in the x-z plane (perpendicular to the Cu02 plane) 
between the Cui3036/Cui2Nd48 and the Cui3 062/Cui2La74 
clusters on a area of 4 a_B x 6 as. The wave functions of 
the cluster electrons have been projected onto the 4s and the 
3d322_^2 states. The solid density lines (black) occur when 
the local electron density of the cluster containing Nd is larger 
than the electron density of the cluster containing La. Dashed 
lines (red), on the other hand, denote that the local electron 
density is larger in the La containing cluster. Two neigh- 
boring lines (solid and dashed) are separated by an electron 
density difference of 0.01 e/sig. 



D. General aspects of Cu EFG 

It is often assumed (by making use of the simple ionic 
model) that the EFG can be used to give an estimate of 
the partial MuUiken population, Pc(3d^2_y2), of the Cu 
3dx^-y2 orbital. This is of course only true if one assumes 
that all the other d orbitals are fully occupied. If the 
population, Pc(3d3z2_,.2), of the Cu 3d^22_j,2 orbital is 
also less than 2, Eq.® must be used for the description of 
the EFG. If we replace the occupation numbers in Eq.|@ 
by the partial MuUiken populations and if we assume 
that ( ) is similar for all 3d orbitals we arrive at the 
following simplified formula: 



C. Comparison with experiments 



(9) 



The nuclear quadrupole frequency, i/q , is related to the 
main component of the EFG, 14^, by 



3e 1 
2hI{2I-l) 



(7) 



Here, Q is the quadrupole moment of the nucleus under 
consideration and rj is the anisotropy parameter. For 
7 = 3/2 and for axial symmetry (ry = 0) one has 



(8) 



Here = pci^d^z^ -r'2) — Pc(3da:2_y2). In Table |IV| 
we have collected the relevant partial MuUiken popu- 
lations Pc{3dx2-y2) and Pci^-^d^z^ ^r^) together with 
and Vzz for several cuprates (including YBa2Cu.30fi 
(Refill), YBazCuaOr (Ref.lH, YBa2Cu408 (Ref.©, 
and Sr2Cu02Cl2 (Ref. For further reference Ta- 

ble IIVI also contains the 4s populations. In Fig. 0] we 
have plotted Vzz versus A^ for various cuprates. The lin- 
ear correlation is evident and we thus conclude that the 
changes in Vzz depend almost entirely on the population 



of the 



and the 3dT,2_j.2 orbitals. This modifi- 



cation of the ionic model tells us that a small value of 
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Vzz can also be obtained by an appropriate population 
difference of the 'iA^z2_^2 and 3d2:2_^2 AOs. 



Substance 


Mulliken populations 


Ad 


v.. 


3d^2_y2 


3d3^2 _^2 


4s 


La2Cu04 


1.406 


1.922 


0.502 


0.516 


1.186 


Nd2Cu04 


1.430 


1.865 


0.647 


0.435 


0.574 


Ndi.g2Ceo.o8Cu04 


1.483 


1.854 


0.644 


0.371 


0.270 


YBaaCuaOr 


1.423 


1.914 


0.590 


0.492 


1.047 


YBa2Cu408 


1.427 


1.909 


0.592 


0.482 


0.986 


YBa2Cu306 


1.428 


1.902 


0.618 


0.474 


0.807 


SraCuOaCla 


1.450 


1.899 


0.525 


0.450 


0.718 



TABLE IV: Partial Mulliken populations of the 3dj.2_y2 and 
3d322_^2 orbitals, their difference Ad, and the calculated EFG 
for La2Cu04, undoped and electron-doped Nd2Cu04, three 
substances of the YBaCuO family, and Sr2Cu02Cl2 for Cui3 
clusters. 



tively. The error bars are estimated from the depen- 
dencies of the size of the cluster. To compare these the- 
oretical values with measurements we need to reassess 
the analysis of orbital shift measurements. As has been 
discussed in detail in Ref. lU the NMR shifts for a nu- 
clear species k are commonly determined as differences of 
frequencies measured in the same magnetic field in the 
target (t) substance and some reference (r) substance. 
Their connection with the chemical shieldings is given 
by 

^Kfit -r) = '^alJr) + ^a^^.Jr) - - 

(10) 

For the copper nucleus, the reference substance is usually 
CuCl and in the above equation the diamagnetic terms 
cancel out: 



63 



Al'(i-CuCl) 



63 „ii 



.(CuCl) 



63 „ii 



At)- (11) 




0.4 0.45 



0.55 



FIG. 4: (color online). EFG versus 

Ad=Pc(3d322_r2)— Pc(3d^2_y2) for La2Cu04 (black square), 
YBa2Cu307 (blue triangle up), YBa2Cu408 (blue triangle 
down), YBa2Cu306 (blue triangle right), Sr2Cu02Cl2 (green 
circle), Nd2Cu04 (red diamond), Ndi.92Ceo.o8Cu04 (full red 
diamond) . 

The relationship between the calculated and the mea- 
sured EFGs will be discussed in depth in a subsequent 
paper—. 



E. Orbital Shifts 

For the cluster Cu50i6/Cu8Nd24 we also have deter- 
mined the chemical shielding tensor for the central cop- 
per nucleus in the same way as has been reported previ- 
ously for La2Cu04 in Ref.|33. The isotropic diamagnetic 
shielding is ^^adia = 2400 ppm. This value is within a few 
ppm the same for all Cu compounds. The paramagnetic 



63, 



para 



(8000 ±800) ppm 
-(2600 ± 300) ppm for the field direc- 
tions parallel and perpendicular to the c-axis, respec- 



shieldings are evaluated as a 
and 



The paramagnetic shielding for CuCl, however, is sub- 
stantial and has experimentally been determined by Lutz 
et al^^ to be ^^a'pZ^iCuClj ^ -1500 ppm with an er- 
ror of about 10 ppm. As it was done in Ref. we in- 
troduce the notation of paramagnetic field modification 

''K''1 = -^uViara- coppcr in Nd2Cu04 we get 



63 



Kl = (0.80 ± 0.08)% and 



s-Vi^ = (0.26 ± 0.03)% 
(12) 

We are not aware of data for orbital shifts in doped 
Nd2Cu04. Detailed measurements, however, have been 
reported by Zheng et alm^ on Pro.giLaCeo.o9Cu04_j,. 

The deduced orbital shifts are ^^kI = 0.92% and 
^^K-^ = 0.155%, measured relatively to the reference 
substance CuCl. After correction for ^^crp^°^(CuCl), as 
discussed above, one obtains the experimental param- 
agnetic field modifications ^'^k\ = 1.07% and ^^AT^ = 
0.305%, respectively. The agreement of the theoretical 
value for the perpendicular component with the experi- 
ment is reasonably good. The predicted value for ^^k\, 
however, is much smaller than that deduced from the 
data. The same discrepancy is also observed for hole 
doped cuprates. We note, however, that the experimental 

value for ^^k\ is based on the observation of a tempera- 
ture independent total magnetic shift and the assumption 
that the on-site and transferred hyperfine fields cancel 
out in Pro.9iLaCeo.o9Cu04_j^, as they do in the YBaCuO 
and LaSCO families. As will be discussed in Sec. lIVI our 
calculated values of the hyperfine couplings give no sup- 
port for these incidental cancellations in substances with 
appreciable differences in the lattice parameters. 



IV. SPIN DENSITY DISTRIBUTION AND 
HYPERFINE COUPLINGS 

The spin-polarized cluster calculations allow us to in- 
vestigate the spin density distribution in detail. In 
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Fig. [S^ we depict the spin density distribution along the 
five Cu and six O obtained for the Cui3 036/Cui2Nd48 
cluster with spin multiplicity M ~ 14. which corresponds 
to a ferromagnetic alignment of the copper moments rep- 
resented in Fig. |5|d. Most of the spin density is provided 
by the singly occupied molecular orbital which is also 
highest in energy. It is a linear combination of idx'^^yi 
atomic orbitals (AOs) on the coppers (~ 80%) and 2po- 
AOs on the oxygens (w 20%). The square of the 'idx^-y2 
AO has maxima at distances of 0.35 A from the nucleus 
but vanishes at the nucleus. In general, the spin density 
at a nuclear site is due to s-type AOs giving rise to the 
Fermi contact term. In particular, the spin densities at 
the Cu nuclei (cusps in Fig.|SJi.) can be correlated with the 
MuUiken spin densities (whose signs are given in Fig.|S|D) 
at the same and the nearest neighbor Cu nuclei. This 
observation can be used to split the Fermi contact term 
into on-site and transferred terms as has been discussed 
in detail in RefsiSSiSS. We just note here that the on- 
site term is negative if the Mulliken spin density at the 
Cu atom under consideration is positive and vice versa 
and that the transferred term is positive (negative) if the 
Mulliken spin densities at the neighboring Cu atoms is 
positive (negative). Near the oxygens, the spin density 
is provided by the 2po- AO and the small cusp at the po- 
sition of the four innermost oxygen nuclei is due to the 
transferred hyperfine fields from the electronic moments 
of the two adjacent copper atoms. In Fig. [St we show 
the spin density distribution obtained with spin multi- 
plicity M = 6. For this multiplicity, the total energy is 
lower than for all other multiplicities and corresponds to 
an antiferromagnetic alignment of the copper moments 
represented in Fig. 

From the differences in the total energy of the various 
spin configurations it is possible to estimate the antifer- 
romagnetic exchange interaction J as will be described 
elsewhera^i. Preliminary results indicate J ~ 160 meV 
for Nd2Cu04 which is about the same as in La2Cu04 
and slightly larger than in YBa2Cu307 (J « 130 meV). 

The hyperfine coupling tensor at the copper is deter- 
mined by an on-site term aj^^ and a transferred term fej^j 
from the nearest neighbors. The value of aj^j is made 
up of three contributions, an isotropic on-site interac- 
tion, aiso, a dipolar interaction, a^'jp, and a contribution 
from spin-orbit coupling, a"^. The transferred term, fo^^j, 
is dominated by the isotropic interaction, hiso, since the 
dipolar term, 6^^*^^, is small and the influence of spin-orbit 
couphng is entirely neglected. In Fig. we show the de- 
pendence of the isotropic hyperfine density, Diso-, at the 
Cu nucleus on the on-site atomic spin density, p^, and 
those of the nearest neighbors, pj. The straight line is a 
fit to the equation 

Aso = aiso • + X] ■ ^^"^^ 

which allows the determination of aiso and hiso- Analo- 
gously, the dipolar contributions can be determined. The 
resulting hyperfine couplings are collected in Table IIVI 



and are compared to the values obtained for La2Cu04. 
The total values a\^^ and include spin-orbit terms 
aio ~ 2.222 a^^ and aj-^ = 0.394 a^"^ which have been es- 
timated according to the procedure given in Ref. We 
note that the values for the two materials are very simi- 
lar. The only significant change is in a^j, which shows a 
different sign in the two substances. 

From experiments there is much less information avail- 
able for the hyperfine couplings in Nd2Cu04 than in the 
hole-doped materials. Zheng et ali^ reported spin shift 
measurements in Pr2_2;Cea;Cu04 and found, with the 
field parallel to the c axis, the same temperature indepen- 
dent behavior as has been observed in doped La2Cu04 
and in the YBaCuO family. The spin shifts measured 
with field in the plane, however, change with decreas- 
ing temperature. This behavior is commonly explained 
by a coincidental cancellation of the hyperfine coupling 
contributions in c direction. The theoretical values for 
aj'^j -I- 46|'qj are indeed close to zero. It is intriguing, 
however, that the same precise cancellation is required 
in most hole doped materials and in Pr2_2;CexCu04 
is intriguing to explain the temperature independent 
spin shifts. More information on hyperfine couplings is 
available for the electron-doped infinite layer compound 
SrCu02. There, -|-46^( should vanish to explain the 
data. Our calculated values give no support for these 
observations. 



ill ' ill ' ill ' ill ' ill ' 

:(a) 1 1 1 1 1 ■ 
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: 1. 1. 1 1 1 1 1 1 1. 1 I : 

l(c) 1 1 1 : 


(d) + 




+ - + 




+ - + - + 
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O Cu O Cu O Cu O Cu O Cu O 



FIG. 5: (color online), (a) Spin density along tfie 5 Cu and 6 
O in the Cui3 036/Cui2Nd48 cluster of spin multiplicity AI = 
14; (c) spin density with multiplicity M = 6. (b) and (d) 
signs of the Mulliken spin densities. 



V. SUMMARY AND CONCLUSIONS 

Spin-polarized ab-initio calculations have been per- 
formed to determine the local electronic structure of 
Nd2Cu04 using clusters comprising 5, 9, and 13 cop- 
per atoms in the Cu02 plane. Electron doping has 
been simulated by two different approaches which both 
yield the same results. The local charge distribution 
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La2Cu04 Nd2Cu04 





-1.94 


-2.42 




-3.55 


-3.38 




0.77 


0.75 




0.08 


0.03 


4ot 


-3.30 


-3.58 


O-tot 


0.22 


-0.34 



TABLE V: Comparative table of the hyperfine constants for 
La2Cu04 and Nd2Cu04 (in units of a^^). 



4 1 1 1 1 1 1 r 




_14l 1 1 1 1 1 1 1 

-8 -6 -4 -2 2 4 6 

r p (Cu.) / p(Cu.) 



FIG. 6: D(Cui)/p(Cui) plotted against Y,j^NN Pi^^j) I 
p(Cui). Values for the cluster with multiplicity M = 14 
(M = 6) are plotted with filled (open) symbols. The differ- 
ent markers correspond to symmetrically inequivalent copper 
sites in the cluster which differ e.g. by their number of nearest 
neighbors. 

has been discussed in detail in terms of the partial Mul- 
liken populations of the individual atomic orbitals. The 



changes in the partial populations of the relevant AOs 
(3d2;2_j,2, 3d322_r2, and 2pcr) that occur upon doping are 
relatively small but significant. Both 'Pci^^x'^-y^^ £^nd 
Pci'ipa) slightly increase and Pc(3d322_^2) is significantly 
less than 2. Electric field gradients depend sensitively on 
the non-spherical local charge distribution. The individ- 
ual contributions to the EFG have been discussed in de- 
tail and the resulting value for the nuclear quadrupole fre- 
quency is in good agreement with the experiments. The 
comparatively small value and the changes upon doping 
have been explained by the differences in the occupation 
numbers N^i^yi and N'^y2_^-2 , which are close to the par- 
tial MuUiken populations Pc- The local electronic struc- 
ture of Nd2Cu04 has been compared to that of La2Cu04. 
Hole doping in the latter material mainly decreases the 
population P(.{idx2-y2) and pd'ipa) but there is also a 
slight decrease of Pc(3d3,2_,.2) and pc(2pz) of the apical 
oxygen. The different values in the EFG at the cop- 
per sites observed in the various cuprates can be under- 
stood by the difference in the partial Mulliken population 
Pc(3d322_r2)— pc(3d^2_j,2). The calculated spin density 
distribution in Nd2Cu04 is not much different from that 
determined for La2Cu04 and YBa2Cu07. The analysis 
of the magnetic hyperfine couplings, which is based on 
the dependence of the contact and dipolar spin densi- 
ties also yields values which are of similar magnitudes to 
those in hole doped materials. 
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